Constraints on 9 13 from A Three-Flavor Oscillation Analysis of 
Reactor Antineutrinos at KamLAND 
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We present new constraints on the neutrino oscillation parameters Am 2 ,!, #12, and #13 from a three-flavor 
analysis of solar and KamLAND data. The KamLAND data set includes data acquired following a radiopurity 
upgrade and amounts to a total exposure of 3.49 x 10 32 target-proton-year. Under the assumption of CPT 
invariance, a two-flavor analysis (#13 = 0) of the KamLAND and solar data yields the best-fit values tan 2 812 = 
0.4441q;°|q and Amli = 7.501q;2o x 10 ~ 5 e V 2 ; a three-flavor analysis with 813 as a free parameter yields 
the best-fit values tan 2 6»i 2 = 0.452j;g;°^, Am 2 ^ = 7.50±g;^ x 1CT 5 eV 2 , and sin 2 8 13 = 0.020+g;^. 
This #13 interval is consistent with other recent work combining the CHOOZ, atmospheric, and long-baseline 
accelerator experiments. We also present a new global #13 analysis, incorporating the CHOOZ, atmospheric 

2 —1-0 013 

and accelerator data, which indicates sin 813 = 0.009_ 007 . A nonzero value is suggested, but only at the 
79%C.L. 

PACS numbers: 14.60.Pq, 26.65,+t, 28.50.Hw, 91.35.-x 
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Neutrino flavor oscillation is by now well established by 
the convergence of results from experiments involving solar, 
reactor, atmospheric, and accelerator neutrinos. Central to 
any discussion of neutrino oscillation phenomenology is the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix 
which describes neutrino mixing in analogy to the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix of the quark sector 11]. 
Although the possibility of more than three neutrino mass 
states, motivated in part by |0], is not excluded, our nota- 
tion and discussion is restricted to the assumption of three 
neutrino mass states. In this case, the three-flavor eigenstates 
(v e , , l> t ) can be expressed as a linear combination of the 
three mass eigenstates {v\ ,vz,v§): 
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Ignoring possible Majorana phases which are irrelevant to os- 
cillation phenomenology, the PMNS matrix U is parametrized 
by three mixing angles, #12, #23, 0i 3 , and a CP-violating 
phase 5. U may be written as 



U 




-iS 



(1) 



rrij) between 



where Sjj = sin0jj and cy — 

The mass-squared splittings (Anijj = m 
the neutrino mass states are described by two independent pa 
rameters, Am^ and Am^. At the currently achieved sensi- 
tivity, mixing between v\ and v-i (M1-V2 mixing) can explain 
the KamLAND data [3] and also, with addition of MSW en- 
hancement iHUt], the solar results Iq-tOl. Atmospheric 11411 . 
K2K and MINOS d] data can be accommodated by 
V2-V3 mixing. As of yet, there is no experimental evidence of 
1/1-1/3 mixing (i.e., a nonzero 6*13) with high statistical signifi- 
cance. 

Probing the value of #13 is a subject of intense ongoing ac- 
tivity. The most stringent limit to date, from the 1 -km-baseline 
CHOOZ reactor experiment (nil, is sin 2 13 < 0.04 at the 
90% C.L. Next-generation accelerator experiments (T2K IU8I1 
and NOj/A [ 19]) and reactor experiments (Double Chooz 12011 . 
Daya Bay 12111 . and RENO i22ll ) aim to significantly improve 
the sensitivity to this parameter and may definitively deter- 
mine the value of #13. If 6*13 is nonzero, future oscillation 
experiments may explore leptonic CP violation (parametrized 
by 8) and probe the neutrino mass hierarchy (i.e., the sign of 
Am| 2 ). The feasibility of such experiments and the path for- 
ward depend critically on the magnitude of #13. 

This article presents an updated KamLAND data set and 
focuses on new constraints on Qyi, Am^, and #13 based on a 
three-flavor combined analysis of KamLAND and solar data. 
As motivated by rf23ll . we also present a global analysis in- 
cluding the CHOOZ, accelerator, and atmospheric oscillation 
experiments in order to explore possible hints of nonzero #13. 



II. APPROXIMATE THREE-FLAVOR NEUTRINO 
OSCILLATION FORMALISM 

Previous KamLAND results were based on a two- 
flavor (y\-V2) oscillation formalism which assumes #13 = 0. 
For the length scales relevant to reactor neutrino oscillation 
at KamLAND and solar neutrino oscillation in the LMA- 
MSW solution, the dependence of the more general three- 
flavor phenomenology on the larger v\-v?, mass splitting 
(|Am 3 i| ~ |Am 32 | 3> Am|x) averages out and the three- 
flavor survival probability (P^), including matter effects, 



may be approximated as 

Pl v = cos 4 
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P e g has the same form as the survival probability in mat- 
ter for V1-V2 mixing but with the electron density (N e ) 
modified: N e — N e cos 2 i3 0. Since sin 2 i3 < 1, 
the survival probability can be further approximated as 
p ee ~ (1 - 2 sin 2 9 13 )?%. Thus, for KamLAND and the 
solar experiments, v\-vz mixing would give rise to an energy- 
independent suppression of the survival probability relative to 
the 6*i3 = case. 

For solar neutrino oscillation in the LMA-MSW solution, 
coherent mixing can be safely ignored due to the long distance 
between the Sun and the Earth. The two-neutrino survival 
probability is simply expressed as 
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where Pf and Pi e are, respectively, the probability of the 
v e — » vi transition in the Sun and the probability of the 
vi — M/ e transition in the Earth with the modified electron den- 
sity N e . Neutrino propagation in the Sun and Earth is calcu- 
lated following the analytical procedure of 11251 12611 . and the 
resulting survival probabilities agree well with numerical cal- 
culations. 

For reactor antineutrinos studied at KamLAND, the matter 
effect in the Earth is not as large as for solar neutrinos. As- 
suming a constant rock density (2.7 g/cm 3 ), the two-neutrino 
survival probability is given by 
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where L is the electron antineutrino (z7 e ) flight distance in me- 
ters from the source to the detector, E is the v e energy in 
MeV, and AmJi is m eV 2 . 0i2M and Am?, 1M are the matter- 
modified mixing angle and mass splitting defined by 
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A = —2\Z2GfN b E, and has a negative sign for antineutri- 
nos; Gf is the Fermi constant. The matter effect modifies the 
expected reactor V e event rate by up to 3%, depending on the 
oscillation parameters. 



III. KAMLAND EXPERIMENT 

The KamLAND detector is located in Kamioka mine, Gifu, 
Japan. The primary target volume consists of 1 kton of ultra- 
pure liquid scintillator (LS). This inner detector (ID) of LS is 
shielded by a 3.2-kton water-Cherenkov outer detector (OD). 
Scintillation light is viewed by 1325 17-inch and 554 20-inch 
photomultiplier tubes (PMTs) providing 34% solid-angle cov- 
erage. A detailed overview of the detector is given in [27]. 
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The T> e flux at KamLAND is dominated by 56 Japanese 
nuclear power reactors. The flux-weighted average baseline 
to these reactors is ^180 km. The reactor fluxes are calcu- 
lated precisely based on detailed operational data including 
the thermal power variation and fuel replacement and reshuf- 
fling records, provided for all Japanese commercial reactors 
by a consortium of Japanese electric power companies. The 
absolute thermal power, used to normalize the fission rates, 
is measured to within 2% for each reactor. This uncertainty 
is conservatively assumed to be correlated across all reactors, 
though some potentially uncorrelated components have been 
put forward in [28]. The data points are typically provided 
at weekly frequency during regular operations when the rel- 
ative instability is of the order of 10~ 3 . When the operating 
parameters vary more quickly, the data are provided at higher 
frequency, with a period between 10 min and 1 h. The rela- 
tive fission yields, averaged over the entire live-time period, 
for isotopes ( 235 U : 238 U : 239 Pu : 241 Pu) are (0.571 : 0.078 
: 0.295 : 0.056), respectively. The detailed reactor operation 
data are also used for accurate tracking of the flux-weighted 
average reactor baseline and spectrum shape change over the 
course of the experiment. The contribution from Korean reac- 
tors, based on reported electric power generation, is estimated 
to be (3.4 ± 0.3)% . The contribution from Japanese research 
reactors and the remainder of the global nuclear power indus- 
try, estimated using reactor specifications from the Interna- 
tional Nuclear Safety Center |2§h is (1.0 ± 0.5)%. The V e 
spectra per fission provided in 030143211 are used, and the un- 
certainties are further constrained from ll33ll . In addition, the 
long-lived, out-of-equilibrium fission products 90 Sr, 106 Ru, 
and 144 Ce [343 are evaluated from the history of fission rates 
for each isotope and are found to contribute only (0.6±0.3)%. 

Electron antineutrinos are detected in KamLAND via the 
inverse beta-decay reaction, V e + p — > e + + n. This process 
has a delayed coincidence (DC) event pair signature which 
offers powerful background suppression. The energy de- 
posited by the positron generates the DC pair's prompt event 
and is approximately related to the incident v e energy by 
E ~ E p + E n + 0.8 MeV, where E p is the sum of the e+ ki- 
netic energy and annihilation 7 energies, and E n is the aver- 
age neutron recoil energy, O(10 keV). The delayed event in 
the DC pair is generated by the capture 7 produced when the 
neutron captures on a proton or 12 C nucleus. The mean neu- 
tron capture time is 207.5 ± 2.8 /is [27] . 

In the previous KamLAND result the largest back- 
ground in the prompt energy region below 3.0 MeV came 
from 13 C(a, n) 16 reactions induced by a-decays in the LS. 
This affected the estimation of the flux of geologically pro- 
duced antineutrinos (geo-F e ) expected between 0.9 MeV and 
2.6 MeV from the decay chains of 238 U and 232 Th in the 
Earth JHUl]. In 2007 the KamLAND collaboration started a 
campaign to purify the LS and ultimately achieved a twenty- 
fold reduction of 210 Po, the dominant a-decay source. This 
reduction gives a better signal-to-background ratio for the 
geo-77 e flux estimation and enhances sensitivity to reactor v e 
oscillations below 2.6 MeV. 

We present an improved measurement of reactor V e oscilla- 
tion based on data collected from March 9, 2002, to November 



4, 2009. This sample includes the previously reported data 
set denoted hereafter as DS-1, in addition to data col- 
lected after LS purification commenced, designated as DS-2. 
The total live time is 2135 days after removing periods of low 
data quality which occurred during LS purification, and af- 
ter detector vetoes to reduce cosmogenic backgrounds. The 
high-quality data selected from DS-2 accounts for 30.41% 
of the total live time. The number of target protons within 
the 6.0-m-radius spherical fiducial volume is calculated to be 
(5.98 ± 0.12) x 10 31 for the combined data set, which corre- 
sponds to an exposure to T> e of 3.49 x 10 32 proton-years. 

Physical quantities such as event vertex and energy are 
reconstructed based on the timing and charge distribu- 
tions of scintillation photons recorded by the ID PMTs. 
The vertex and energy reconstructions are calibrated us- 
ing 60 Co, 68 Ge, 203 Hg, 65 Zn, 241 Am 9 Be, 137 Cs, and 
2ln Po 13 C radioactive sources. The observed vertex reso- 
lution is ~12 cm/ y/ E(MeV), and the energy resolution 
is 6.4%/\/i?(McV). For DS-2, the resolutions are time- 
dependent due to a light-yield reduction of up to ^20% rel- 
ative to DS-1. The source calibrations are augmented with 
studies of muon spallation products to monitor the detector 
stability and to determine the nonlinearity of the energy re- 
sponse due to LS quenching, Cherenkov light, and dark hit 
contributions. The systematic uncertainty of the absolute en- 
ergy response over the full DS-1 and DS-2 data sets is less 
than 1.2%, and when propagated in the reactor V e spectrum 
produces a 1.8% uncertainty on Am^ and a 1.3% uncertainty 
on the event rate above the analysis threshold. 

For DS-1, the systematic uncertainty on the fiducial vol- 
ume up to 5.5 m radius was determined to be 1.6% with a 
full-volume calibration campaign 13711 . The uncertainty in 
the volume between 5.5 m and 6.0 m radius was estimated 
from the vertex uniformity of muon-induced 12 B and 12 N; the 
combined uncertainty on the 6.0-m-radius fiducial volume for 
DS-1 is 1.8%. To date there have been no full-volume calibra- 
tions for DS-2, so we rely on vertex uniformity of cosmogenic 
12 B and 12 N events; in this case, we assign a 2.5% uncertainty 
on the 6.0-m-radius fiducial volume. 

Table Q] summarizes the systematic uncertainties on Am^ 
and the expected event rate of reactor 77 e 's; the overall rate 
uncertainties for DS-1 and DS-2 are 4.1% and 4.5%, respec- 
tively. 



IV. KAMLAND DATA REDUCTION AND CANDIDATE 
EVENT SELECTION 

Antineutrino DC-pair candidates are selected by perform- 
ing the following series of first-level cuts: (i) prompt en- 
ergy: 0.9 < £p(MeV) < 8.5 ; (ii) delayed energy: 1.8 < 
E d (MeV) < 2.6 (capture on p), or 4.4 < E d (MeV) < 5.6 
(capture on 12 C); (iii) spatial correlation of prompt and de- 
layed events: Ai?(m) < 2.0 ; (iv) time separation between 
prompt and delayed events: 0.5 < AT(/is) < 1000 ; and (v) 
fiducial volume radii: i? p ,i?d(m) < 6.0 . 

In order to increase the ratio of signal to accidental- 
background, a second-level cut is performed using a likeli- 
hood discriminator, L = ,_ J "° — . Here f Ve and f acc are the 
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TABLE I: Estimated systematic uncertainties for the neutrino oscil- 
lation parameters Am|i, #12, and #13 for the earlier/later periods of 
measurement, denoted in the text as DS-l/DS-2. The overall uncer- 
tainties are4.1% / 4.5% for DS-l/DS-2. 



Detector-related (%) 



Reactor-related (%) 



Rate 



Energy scale 

Fiducial volume 
Energy scale 
Lcut(E p ) eff. 
Cross section 
Total 



1.8/1.8 

1.8/2.5 
1.1/1.3 
0.7/0.8 
0.2/0.2 
2.3/3.0 



^ e -spectra [33] 

77 e -spectra 
Reactor power 
Fuel composition 
Long-lived nuclei 
Total 



0.6/0.6 
2.4/2.4 
2.1/2.1 
1.0/1.0 
0.3/0.4 
3.3/3.4 



probability density functions (PDFs) for v e DC pairs and ac- 
cidental DC pairs, respectively; both PDFs are functions of 
the 6 DC-pair parameters: E p , E d , AR, AT, R p , R d . The 
PDF for accidental DC pairs can be evaluated directly from 
the data with an off-time cut; we use 10 ms < AT < 20 s. To 
utilize the variation in the accidental DC rate with time, the 
full data set is divided into five periods and the corresponding 
f acc is computed for each. The PDF for 77 e DC pairs is cal- 
culated with a Monte Carlo (MC) simulation. The systematic 
error in the simulated PDF is evaluated by comparing simu- 
lated calibration data to real calibration data for the 68 Ge and 
241 Am 9 Be sources. 

For each 0.1 MeV interval in prompt energy, we choose 
L C ut(E p ) to maximize ^ S + Ba 



where S and B acc are the 
expected number of 77 e and accidental DC pairs, respectively, 
with L(E p ) > L cut (E p ). To exploit the time variation of 
both the signal and background, the optimal L cut (E p ) is de- 
termined for each of the five time periods. Finally, only DC 
pairs with L(E p ) > L cut (E p ) are selected. The efficiency 
and uncertainty of the cut are evaluated for each period us- 
ing the MC; the independent efficiency, averaged over the 
five time periods, is shown in the top panel of Fig. Q] A no- 
oscillation input spectrum is used to generate . The effect 
of using an oscillated v e spectrum was checked with various 
trial values of (#12, Ato 2 i) an d found not to greatly affect the 
selection. The number of accidental DC pairs remaining af- 
ter all cuts is determined to be 102.5 ± 0.1. The dominant 
contributors to these accidental DC pairs are 2.6 MeV 7-rays 



TABLE II: Estimated backgrounds excluding geo-^ e after first- and 
second-level cuts. 



Background 



Contribution 



Accidental 
9 Li/ 8 He 
13 C(a,n) 16 
13 C(a,n) 16 O g . s ., rA C(n,n 
13 C(q, n) ie O, 1st e.s. (6.05 MeV e + e") 
13 C(q, n) ie O, 2nd e.s. (6.13 MeV 7) 
Fast neutron and atmospheric neutrino 



np — > np 

i2rv„ „'U2q» (44 MeV 7) 



Total 



102.5 ± 0.1 

24.8 ± 1.6 
171.7 ± 18.2 

7.3 ± 0.8 

15.9 ± 3.3 
3.7 ± 0.7 
< 12.3 

325.9 ± 26.1 
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FIG. 1 : Prompt energy spectrum of 77 e candidate events above 0.9 
MeV energy threshold (vertical dashed line). The data together with 
the background and reactor z7 e contributions fitted from an unbinned 
maximum-likelihood three-flavor oscillation analysis are shown in 
the main panel. The number of geo-F e 's is unconstrained in the fit. 
The shaded background histograms are cumulative. The top panel 
shows the energy-dependent selection efficiency; each point is the 
weighted average over the five time periods described in the text. 



from external 208 T1 /3-decays. 

In addition to accidental background events, there are 
other processes which produce background DC pairs. The 
13 C(a, n) 16 nuclear reaction in the LS is the largest such 
background. The dominant a source is 210 Po, a long-lived 
daughter nucleus of 222 Rn. This reaction produces neutrons 
with energies up to 7.3 MeV, and mostly contributes DC 
pairs with prompt energies below 2.6 MeV. By counting the 
quenched scintillation signals from the 5.3 MeV a particles, 
we find (5.95 ± 0.29) x 10 9 a-decays in full data set. The 
rate of the 13 C(a, n) 16 background and its prompt energy 
spectrum is estimated by simulation. The total cross sec- 
tion and final-state partial cross sections for 16 O, Oi (where 
i = 0,1,2 for the ground, first and second excited states 
of 16 0), are based on |[38l [3911 . but the relative normaliza- 
tions of the <Ji were tuned by an in-situ calibration using a 
210 Po 13 C source 14011 . The data require ctq and o\ be scaled 
by 1.05 and 0.6, respectively, while no scaling is required for 
(72 • Including the uncertainty on the number of a-decays, we 
assign an uncertainty of 1 1 % for the ground state and 20% 
for the excited states. We estimate that the total number of 
13 C(a, n) 16 DC pairs remaining in the full data set after 
the first- and second-level cuts is 198.6 ± 23.0. DS-2, which 
benefited from reduced 210 Po contamination due to LS purifi- 
cation, contributes only 7% of the 13 C(a, n) 16 events after 
all selection cuts. 
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FIG. 2: Allowed regions projected in the (tan 2 612, Amfx) plane, for solar and KamLAND data from (a) the two-flavor oscillation analysis 
(#13 = 0) and (b) the three-flavor oscillation analysis, where 813 is a free parameter. The shaded regions are from the combined analysis of 
the solar and KamLAND data. The side panels show the A\ 2 profiles projected onto the tan 2 #12 and Am^ axes. 



Delayed-neutron beta emitters 9 Li and 8 He, which are pro- 
duced in the LS by cosmic-ray muons, also generate DC 
pairs 112711 . They are removed by a 2-s veto of the entire fidu- 
cial volume after LS showering muons, which generate more 
than 10 6 photoelectrons in the LS, and poorly reconstructed 
LS muons. In the case of nonshowering, well-reconstructed 
LS muons, the 2-s veto is applied only within a 3-m-radius 
cylinder around the muon track in order to minimize the ex- 
posure loss from the veto. From a fit to the time delay between 
prompt DC events and their preceding LS muons, we estimate 
the background remaining after the veto and DC selection cuts 
is 24.8 ± 1.6 events. 

Fast neutrons and atmospheric neutrinos are also a possible 
source of DC pairs. Fast neutrons generated in the material 
outside the OD may scatter into the ID, and subsequent co- 
incidence signals in the LS from prompt neutron scatter and 
delayed capture sometimes pass the v e DC signal selection 
criteria. Monte Carlo studies of neutron generation outside 
the ID 1 27] indicate that fast neutrons are generated primarily 
by cosmic-ray muons. A 2-ms veto after OD-tagged muons 
mostly eliminates fast neutron DC pairs. The residual back- 
ground due to the OD tagging inefficiency and muons that 
pass nearby but do not enter the OD is estimated from simu- 
lation. Atmospheric neutrino backgrounds are evaluated us- 
ing the NUANCE software 114 ill to simulate neutrino interac- 
tions and related processes. Both atmospheric neutrino and 
fast neutron DC pairs are assumed to have a flat prompt en- 
ergy spectrum in the energy range of the present analysis, and 
are estimated to contribute less than 12.3 candidates in total 
after all selection cuts. 

Geo-T7 e fluxes at Kamioka can be calculated based on a ref- 
erence Earth model 04211 which assumes a radiogenic heat pro- 
duction rate of 16 TW from the decay chains of U and Th. 



Including neutrino oscillation effects, this model predicts 85 
and 21 events in the full data set from U and Th decays, re- 
spectively. However, since the estimate of the geo-77 e yield is 
highly dependent on the Earth model, the event rates from the 
U and Th decay chains are not constrained in the oscillation 
analysis; only the prompt energy spectrum shapes, which are 
independent of the Earth model, are used to constrain their 
contribution. A possible background contribution from a hy- 
pothetical reactor-I7 e source at the Earth's center, motivated 
by 114311 and investigated in H36I1 . is neglected in this analysis. 

After all selection cuts, we expect, in the absence of v e 
disappearance, 2879±118 events from reactor V e , and 325. 9± 
26.1 events from the backgrounds, as summarized in Table [Hi 
The observed number is 2106 events. 



V. OSCILLATION ANALYSIS 

The KamLAND data is analyzed based on an unbinned 
maximum-likelihood method. The \ 2 is defined by 

X 2 = X?atc( 6 'l2,^13, Amli, NBG1^5, Nf^Th'"!-^) 

-21nL shap c(6 , i2, #13, Am^, A^bgi^s, N^ h , cki-ml) 

+XBG (^BG1^5) + XsystOl^O ■ (7) 

The terms are, in order: the \ 2 contribution for (i) the to- 
tal rate, (ii) the prompt energy spectrum shape, (iii) a penalty 
term for backgrounds, and (iv) a penalty term for systematic 
uncertainties. A^bgi->5 are tne expected background levels 
discussed in Sec. [IVj and N^ c< ^ h are the contributions ex- 
pected from U and Th geo-I7 e 's. -/Vbgi->5 Slle allowed to 
vary in the fit but are constrained with the penalty term (iii) 
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FIG. 3: Allowed regions from the solar and KamLAND data pro- 
jected in the (tan 2 #12, sin 2 #13) plane for the three-flavor analysis. 




FIG. 4: Ax 2 -profiles projected onto the sin 2 #13 axis for different 
combinations of the oscillation data floating the undisplayed param- 
eters (tan 2 612, Am^i). 



using the estimates summarized in Table HH -ZVy™ h are free 
parameters and are unconstrained to avoid any Earth model 
dependence. The 01^4 parametrize the uncertainties on the 
reactor 77 e spectra and energy scale, the event rate, and the 
energy dependent efficiencies; these parameters are allowed 
to vary in the analysis but are constrained by term (iv). The 
background energy scale uncertainties are estimated to con- 
tribute at most an additional 0.5% to the error on the event 
rate and are neglected in this analysis. The prompt energy 
spectrum shape likelihood term (ii) is evaluated as a func- 
tion of the candidate event time. The detailed knowledge 
of the time evolution of the total reactor v e spectrum and 
effective baseline, afforded by the reactor fuel composition 
and power data provided by the Japanese reactor operators, 
is thus fully utilized in the analysis. Variations in the to- 
tal observed spectrum shape with time due to changes in 
the background levels — especially the 13 C(a, n) 16 reduc- 
tion from the LS purification — are also exploited by this term. 
The spectrum shape likelihood term allows an Earth-model- 
independent constraint of the geo-77 e contribution since the U 
and Th decay spectra are known independently of the Earth 
model. A x 2 scan of the (#12, #13, Am 2 ^) oscillation pa- 
rameter space is carried out, minimizing x 2 with respect to 
N B gi->5, ^u™h' and 0^4- 

In our analysis of the solar neutrino data, we include 
the rates in the chlorine [6] and gallium [9] experiments, 
Borexino d], SNO III QJ, the zenith spectra in Super- 
Kamiokande phase I iflfjll . and the day-night spectra in SNO 
phase I and II Hfl]. 

The measured fluxes are compared 
with the high-metallicity standard solar model predictions 
(GS98) H. 

For the three-flavor KamLAND-only analysis, without 
any constraints on #13 from other oscillation experiments, 



the best-fit oscillation parameter values are Am^ — 
7.49±g;|g x 10~ 5 eV 2 , tan 2 8 12 = 0.436±[J;^ 2 and 
sin 2 #13 = 0.0321°;°^ « 0.094 at the 90% C.L.). The two- 
flavor oscillation treatment using Eq. (|7]), as presented pre- 
viously in [3], is a special case of the three-flavor treat- 
ment with 6*i3 = 0. For this case the best-fit oscillation 
parameters from the KamLAND-only analysis are Am 2 ^ = 
7.50+°;2o x 10" 5 eV 2 and tan 2 6 12 = 0.492+°;°^. In the 
KamLAND data, #13 is expected to contribute only an energy- 
independent event rate suppression and we find almost no ef- 
fect on the Ato^ measurement when #13 is included as a 
free parameter. Figure [T] shows the prompt energy spectrum 
of candidate events in KamLAND together with the best-fit 
background and reactor 77 e spectra for the three-flavor fit to 
the KamLAND data. The fit estimates 82 and 26 events from 
U and Th geo-77 e 's, respectively, in agreement with the refer- 
ence model. 

Figure |2] compares the allowed regions in the 
(tan 2 #12, Arn 2 i) plane from the two- and three-flavor 
oscillation analyses. We find [Fig. Eta)] that the al- 
lowed region from the solar data is in agreement with the 
KamLAND data, and the small tension between the two- 
flavor best-fit values of 812, discussed previously in l23ll45ll . 
has eased. Assuming CPT invariance, the two-neutrino 
oscillation parameter values from a combined analysis of 
the solar and KamLAND data are tan 2 8 12 = 0.444±o;o3o 
and Am 2 .! = 7.50+^0 x 10~ 5 eV 2 . For the three-flavor 
analysis combining the solar and KamLAND data, the 
best-fit parameter values are tan 2 812 = 0.452^0 033 and 
sin 2 013 = 0.020l^^; the best-fit value for Am| x is the 
same as for the two-flavor result. The best-fit values for 
the different data combinations and analysis approaches are 
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summarized in Table HHI in Appendix A. 

Figure[3]shows the regions in the (tan 2 # 12 , sin 2 0i 3 ) plane 
allowed by \ 2 minimization with respect to Am|i for each 
analysis. The reduction of the best-fit value of tan 2 0i2 for 
the three-flavor KamLAND-only analysis relative to the two- 
flavor KamLAND analysis (Fig. |2} follows the anticorre- 
lation apparent in the KamLAND contours (Fig. [3). The 
correlation between 6*12 and #13 in the solar data is slight and 
the difference between the best-fit values of 612 from the two- 
flavor and three-flavor analyses of the solar-only data is small. 

Figure |4] shows Ax 2 -profiles projected onto the sin 2 #13 
axis for different combinations of the data. The analysis of 
the KamLAND data gives sin 2 8 13 = 0.032±g;g|^ (< 0.094 
at the 90% C.L.), and the combined analysis of the solar and 
KamLAND data gives sin 2 6 13 = 0.020ig;[£{j. The con- 
straint on nonzero #13 from the combined KamLAND and 
solar analysis is comparable to the constraint from the com- 
bined analysis of CHOOZ, atmospheric, and long-baseline ac- 
celerator (LBL, i.e., K2K and MINOS) experiments presented 
in l46ll . which includes the recent v e appearance result from 
MINOS H. In the solar + KamLAND analysis the pref- 
erence for nonzero #13 comes mostly from the KamLAND 
data. All oscillation data favor a positive 6*13, although the 
current statistical power is poor. For a global analysis com- 
bining our updated KamLAND + solar analysis with the com- 



bined CHOOZ, atmospheric, and LBL (appearance + disap- 

^^^^t^ o _i_o ni q 

pearance) analysis from 114611 . we find sin 6*13 = Q.009_q 007. 
Our global result is very similar to the global analysis carried 
out by [46] with the previous KamLAND data set but the sig- 
nificance of nonzero #13 is reduced slightly to the 79% CL. 



VI. VISUALIZATION OF THE SURVIVAL PROBABILITY 

Figures [5]and |6]illustrate different aspects of the survival 
probability for the KamLAND data. The data points in Fig. 
[5] are the ratio of the observed reactor v e spectrum to that ex- 
pected in the case of no oscillation plotted as a function of 
Lq/E, where Lq (Lq = 180 km) is the flux-weighted average 
reactor baseline. The oscillatory structure arising from the 
sin 2 (1.27 Am^L/E) term is clear, but is distorted because 
the reactor sources are distributed across multiple baselines. 
We also overlay in the figure the expected oscillation curves 
based on the best-fit parameters from the two- and three-flavor 
unbinned maximum-likelihood analyses discussed previously. 
The suppression of the oscillation amplitude is slightly larger 
for the nonzero 6*13 case. 

To focus on 6*12 and #13 effects in the data, we introduce a 
parameter x(E p , t) defined by 



x(E p ,t) 



sin 2 20i2 



N n 



E 



dE sin 2 2f?i2 M sin 2 f L27A ™2iM^ 



E 



1 



sin 2 20i2 



.2™ ■ 2 ( L27AmS 1M L 
sin 2 2012M sin 2 f g 21M 



(8) 
(9) 



where 

reactors „ Q ( T? +\ 

N no . osc (E p ,t) = I dEP *(E P ,t,E)^-^ (10) 

is the number of candidates with prompt energy E p expected 
in the absence of neutrino oscillation from all reactors at time 
t at KamLAND; the index i labels the reactor source; Li and 
Si(E, t) are, respectively, the baseline and the neutrino spec- 
trum at time t of reactor i; and Pr(E p , t, E) is the probability 
that a V e with energy E will be detected at KamLAND with 
prompt energy E p . Pr includes the number of target protons, 
the inverse beta-decay cross section, and the time-dependent 
detector response function. ($12, Am^) are the best-fit values 
from the two-flavor unbinned analysis, and (0i2A/,Am| 1Af ) 
are the matter-modified oscillation parameters calculated with 
those best-fit values. The angle bracket notation in Eq. (0 
indicates the weighted average over reactor baselines Li and 
neutrino emission energies E, written explicitly in Eq. (0. 
For the region of (Am|i, #12, $13) parameter space close to 
Am|i, all the information about the reactors, detector-related 
effects, and matter modification is contained in the parameter 



x. With this definition, the survival probability may be written 
as a linear function of x, P(E p , t) = A — B ■ x(E p ,t) , where 
A = (cos 4 6»i3 + sin 4 13 ) and B = cos 4 i3 sin 2 29 12 . 9 13 
effects are predominantly encoded in A, whereas 6*12 effects 
dominate the slope B. This linear relationship is illustrated 
in Fig. [6] The points there are the survival probability for 
KamLAND events binned as a function of x. Also shown 
are lines where A and B have been calculated using the best- 
fit values from the two- and three-flavor unbinned maximum- 
likelihood analyses of the KamLAND data. The axis intercept 
at x = of the best-fit 3-v line is less than one, illustrating the 
slight indication of positive 6*13 from the unbinned likelihood 
analysis. Any further improvement in the significance of the 
013 investigation with KamLAND requires reduced system- 
atic uncertainties on the reactor flux and increased detector 
exposure. A binned analysis based on the data points in Fig. 
[6] is outlined in Appendix B. 
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FIG. 5: Ratio of the observed v e spectrum to the expectation for no- 
oscillation versus Lo/E for the KamLAND data. Lq = 180 km 
is the flux-weighted average reactor baseline. The 2-v and 3-v his- 
tograms are the expected distributions based on the best-fit parameter 
values from the two- and three-flavor unbinned maximum-likelihood 
analyses of the KamLAND data. 



VII. CONCLUSION 

An updated KamLAND reactor v e data set was presented. 
The data set benefits from increased exposure and an im- 
proved background environment due to a radiopurity upgrade 
of the LS. The analysis slightly hints at a nonzero 6\% with the 
available oscillation data. In a two-flavor analysis (0 13 = 0) 
of the solar and KamLAND data, the best-fit values for the os- 
cillation parameters are tan 2 6*12 = 0.444to'o3o an ^ ^ m 2i = 
7.50^20 x 10~ 5 eV 2 . In the three-flavor analysis, floating 
the value of $13 without any constraints from the other oscil- 
lation experiments gives the solar + KamLAND best-fit values 
tan 2 0i2 = 0.452i°:°33, Amfi = 7.50±£ ^ x 1(T 5 eV 2 , and 
sin 2 013 = 0.020lQQig. The limits on Amfi are the same for 
the two- and three-flavor analyses. All three oscillation pa- 
rameters derived from the KamLAND-only antineutrino data 
are in good agreement with those from the solar-only neu- 
trino data and reveal no inconsistency with CPT invariance, 
which was assumed for the joint fits. The upper limit we ob- 
tain on sin 2 #13 is compatible with other recent work combin- 
ing CHOOZ, atmospheric, and accelerator experiments. More 
definitive information on the value of 6*13 should come from 
upcoming accelerator and reactor experiments. 
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APPENDIX A 

The best-fit values for the different data combinations and 
analysis approaches are summarized in Table ITTT1 



APPENDIX B 

We consider the unbinned maximum-likelihood method 
presented in Sec. [V] to be the optimal approach to analyz- 
ing the KamLAND data because it takes full advantage of all 
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TABLE III: Summary of the best-fit values for tan 2 6*12 and sin 2 #13 
from two- and three-flavor neutrino oscillation analyses of various 
combinations of experimental data. 



Data set 



Analysis method tan #12 



KamLAND 
KamLAND + solar 
KamLAND 
KamLAND + solar 
Global 



two-flavor 
two- flavor 
three-flavor 
three-flavor 
three-flavor 



400+0 086 
444 +o.036 

0.436 



0.030 
+0.102 



-0.081 

u '™-0 .033 
+0.039 
0.027 



0.444 



= 
= 

U.UOZ_g q 37 

o.o2o±8:gi| 

0.0091S-SS 



the spectral and time information available. In this appendix 
we outline a binned oscillation analysis which we find repro- 
duces very well the A\ 2 contours in the (6*12, #13) subspace 
for the unbinned likelihood KamLAND-only analysis shown 
in Fig. [3] The binning parameter is the parameter x intro- 
duced in Sec. |VT] and defined in Eq. ®. Table [TV] lists the 
binned data. The binned \ 2 is defined as 



E 



Pi - Pi(l + 5 corT ) 



where 



1 



sin 2 26*12 ■Xi) + sin 4 6>i 3 



, (ID 



(12) 



and Xi is the weighted average of x over bin i. The pairs 
(pi, a Pi ) are the observed survival probability, defined as the 
ratio of the observed events to the expectation for no os- 
cillation, and its uncertainty for each bin i, and 6 COTT is 
a factor needed to account for the systematic uncertainty 



(fcorr = 0.041) on the flux prediction. In Eq. ( fl~2b . the vac- 
uum 8 12 should be used because matter corrections to 812 and 
Am 2 ! are included in the calculation of x, as shown in Eq. ©. 
For a global analysis, the small dependence on Am^ can be 
ignored and the binned \ 2 mav be use d f° r a scan over the 
(012,013) oscillation parameter space. Comparing the Ax 2 
map built using this method and that from the full unbinned 
analysis shown in Fig. [3] the only significant deviations ap- 
pear far from the best-fit point at high values of 812 where 
constraints from the solar neutrino experiments dominate. 

TABLE IV: Survival probability for each bin in x [de- 
fined in Eq. {8}]. The first column indicates the bin range 
of ie (sin 2 2012m sin 2 (1.27 Aml 1M L/E)) / sin 2 2§ 12 . The 
weighted average x is given in the second column. The observed 
survival probability is shown in the third column. The uncertainties 
include only the statistical and background estimation uncertainties, 
which are assumed to be uncorrelated. In addition, the systematic un- 
certainty (cr corr = 4.1%) on the flux prediction needs to be included 
for each bin as a fully correlated uncertainty. 



x range 


X 


Survival probability (p ± a p ) 


0.1-0.3 


0.230 


0.749 ± 0.044 


0.3-0.4 


0.354 


0.650 ± 0.039 


0.4-0.5 


0.451 


0.624 ± 0.046 


0.5-0.6 


0.555 


0.512 ±0.038 


0.6-0.7 


0.638 


0.416 ± 0.030 


0.7-0.9 


0.800 


0.415 ±0.160 



[1] K. Nakamura et al. (Particle Data Group), J. Phys. G 37, 075021 
(2010). 

[2] A. Aguilar, L. B. Auerbach, R. L. Burman, D. O. Caldwell, 
E. D. Church, A. K. Cochran, J. B. Donahue, A. Fazely, G. T. 
Garvey, R. M. Gunasingha, et al. (LSND Collaboration), Phys. 
Rev. D 64, 112007 (2001). 

[3] S. Abe, T. Ebihara, S. Enomoto, K. Furuno, Y. Gando, 
K. Ichimura, H. Ikeda, K. Inoue, Y. Kibe, Y. Kishimoto, 
et al. (KamLAND Collaboration), Phys. Rev. Lett. 100, 221803 
(2008). 

[4] L. Wolfenstein, Phys. Rev. D 20, 2634 (1979). 
[5] S. P. Mikheyev and A. Y. Smirnov, Sov. J. Nucl. Phys. 42, 913 
(1985). 

[6] B. T. Cleveland, T. Daily, R. Davis, Ir., I. R. Distel, K. Lande, 
C. K. Lee, P. S. Wildenhain, and J. Ullman, Astrophys. J. 496, 
505 (1998). 

[7] W. Hampel, I. Handt, G. Heusser, I. Kiko, T. Kirsten, 

M. Laubenstein, E. Pernicka, W. Rau, M. Wojcik, Y. Zakharov, 

et al., Phys. Lett. B 447, 127 (1999). 
[8] M. Altmann, M. Balata, P. Belli, E. Bellotti, R. Bernabei, 

E. Burkert, C. Cattadori, R. Cerulli, M. Chiarini, M. Cribier, 

et al., Phys. Lett. B 616, 174 (2005). 
[9] I. N. Abdurashitov, V. N. Gavrin, V. V. Gorbachev, P. P. Gurk- 

ina, T. V. Ibragimova, A. V. Kalikhov, N. G. Khairnasov, T. V. 

Knodel, I. N. Mirmov, A. A. Shikhin, et al., Phys. Rev. C 80, 



015807 (2009). 

[10] I. Hosaka, K. Ishihara, I. Kameda, Y. Koshio, A. Minamino, 
C. Mitsuda, M. Miura, S. Moriyama, M. Nakahata, T. Namba, 
et al. (Super-Kamiokande Collaboration), Phys. Rev. D 73, 
112001 (2006). 

[11] B. Aharmim, S. N. Ahmed, A. E. Anthony, N. Barros, E. W. 
Beier, A. Bellerive, B. Beltran, M. Bergevin, S. D. Biller, 
K. Boudjemline, et al. (SNO Collaboration), Phys. Rev. C 81, 
055504 (2010). 

[12] B. Aharmim, S. N. Ahmed, I. F. Amsbaugh, A. E. Anthony, 
I. Banar, N. Barros, E. W. Beier, A. Bellerive, B. Beltran, 
M. Bergevin, et al. (SNO Collaboration), Phys. Rev. Lett. 101, 
111301 (2008). 

[13] C. Arpesella, H. O. Back, M. Balata, G. Bellini, I. Benziger, 
S. Bonetti, A. Brigatti, B. Caccianiga, L. Cadonati, F. Calaprice, 
et al. (Borexino Collaboration), Phys. Rev. Lett. 101, 091302 
(2008). 

[14] R. Wendell, C. Ishihara, K. Abe, Y. Hayato, T. Iida, M. Ikeda, 
K. Iyogi, I. Kameda, K. Kobayashi, Y. Koshio, et al. (Super- 
Kamiokande Collaboration), Phys. Rev. D 81, 092004 (2010). 

[15] M. H. Ahn, E. Aliu, S. Andringa, S. Aoki, Y. Aoyama, I. Ar- 
gyriades, K. Asakura, R. Ashie, F. Berghaus, H. G. Berns, et al. 
(K2K Collaboration), Phys. Rev. D 74, 072003 (2006). 

[16] P. Adamson, C. Andreopoulos, K. E. Arms, R. Armstrong, D. I. 
Auty, D. S. Ayres, B. Bailer, P. D. Barnes, G. Barr, W. L. 



10 



Barrett, et al. (MINOS Collaboration), Phys. Rev. Lett. 101, 
131802 (2008). 

[17] M. Apollonio, A. Baldini, C. Bemporad, E. Caffau, F. Cei, 

Y. Declais, H. d. Kerret, B. Dieterle, A. Etenko, L. Foresti, et al., 

Eur. Phys. J. C 27, 331 (2003). 
[18] Y. Itow, T. Kajita, K. Kaneyuki, M. Shiozawa, Y. Totsuka, 

Y. Hayato, T. Ishida, T. Ishii, T. Kobayashi, T. Maruyama, et al. 

(2001), hep-ex/0106019vl. 
[19] D. Ayres et al. (2005), hep-ex/0503053vl. 
[20] F. Ardellier, I. Barabanov, J. C.Barriere, M. Bauer, L. Bezrukov, 

C. Buck, C. Cattadori, B. Courty, M. Cribier, F. Dalnoki-Veress, 

etal. (2004), hep-ex/0405032vl. 
[21] X. Guo et al. (2007), hep-ex/0701029vl. 
[22] J. K. Ahn et al. (2010), hep-ex/1003.1391vl. 
[23] G. L. Fogli, E. Lisi, A. Marrone, A. Palazzo, and A. M. Ro- 

tunno, Phys. Rev. Lett. 101, 141801 (2008). 
[24] S. Goswami and A. Y. Smirnov, Phys. Rev. D 72, 053011 

(2005). 

[25] S. J. Parke, Phys. Rev. Lett. 57, 1275 (1986). 

[26] E. Lisi and D. Montanino, Phys. Rev. D 56, 1792 (1997). 

[27] S. Abe, S. Enomoto, K. Furuno, Y. Gando, H. Ikeda, K. In- 

oue, Y. Kibe, Y. Kishimoto, M. Koga, Y. Minekawa, et al. 

(KamLAND Collaboration), Phys. Rev. C 81, 025807 (2010). 
[28] Z. Djurcic, J. A. Detwiler, A. Piepke, V. R. Foster, L. Miller, 

and G. Gratta, J. of Phys. G: Nucl. and Part. Phys. 36, 045002 

(2009). 

[29] INSC, the International Nuclear Safety Center available at 

http://www.insc.anl.gov/ (2010). 
[30] K. Schreckenbach, G. Colvin, W. Gelletly, and F. Von Feil- 

itzsch, Phys. Lett. B 160, 325 (1985). 
[31] A. A. Hahn, K. Schreckenbach, W. Gelletly, F. von Feilitzsch, 

G. Colvin, and B. Krusche, Phys. Lett. B 218, 365 (1989). 
[32] P. Vogel, G. K. Schenter, F. M. Mann, and R. E. Schenter, Phys. 

Rev. C 24, 1543 (1981). 
[33] B. Achkar, R. Aleksan, M. Avenier, G. Bagieu, J. Bouchez, 



R. Brissot, J. F. Cavaignac, J. Collot, M.-C. Cousinou, J. P. Cus- 

sonneau, et al., Phys. Lett. B 374, 243 (1996). 
[34] V. I. Kopeikin, L. A. Mikaelyan, and V. V. Sinev, Phys. of At. 

Nucl. 64, 849 (2001). 
[35] T. Araki, S. Enomoto, K. Furuno, Y. Gando, K. Ichimura, 

H. Ikeda, K. Inoue, Y. Kishimoto, M. Koga, Y. Koseki, et al. 

(KamLAND Collaboration), Nature 436, 499 (2005). 
[36] G. Bellini, J. Benziger, S. Bonetti, M. B. Avanzini, B. Cac- 

cianiga, L. Cadonati, F. Calaprice, C. Carraro, A. Chavarria, 

F. Dalnoki-Veress, et al. (Borexino Collaboration), Phys. Lett. 

B 687, 299 (2010). 
[37] B. E. Berger, J. Busenitz, T. Classen, M. P. Decowski, D. A. 

Dwyer, G. Elor, A. Frank, S. J. Freedman, B. K. Fujikawa, 

M. Galloway, et al., J. Instrum. 4, P04017 (2009). 
[38] S. Harissopulos, H. W. Becker, J. W. Hammer, A. Lagoyannis, 

C. Rolfs, and F. Strieder, Phys. Rev. C 72, 062801 (2005). 
[39] JENDL, the Japanese Evaluated Nuclear Data Library available 

at http://wwwndc.jaea.go.jp/jendl/jendl.html (2005). 
[40] D. W. McKee, J. K. Busenitz, and I. Ostrovskiy, Nucl. Instrum. 

and Meth. A 587, 272 (2008). 
[41] D. Casper, Nucl. Phys. B - Proc. Supp. 112, 161 (2002). 
[42] S. Enomoto, E. Ohtani, K. Inoue, and A. Suzuki, Earth and 

Planet. Sci. Lett. 258, 147 (2007). 
[43] J. M. Herndon, Proc. Natl. Acad. Sci. U.S.A. 100, 3047 (2003). 
[44] N. Grevesse and A. J. Sauval, Space Sci. Rev. 85, 161 (1998). 
[45] M. Gonzalez-Garcia, M. Maltoni, and J. Salvado, J. High En- 
ergy Phys. 4, 056 (2010). 
[46] M. Gonzalez-Garcia, M. Maltoni, and J. Salvado (2010), hep- 

ph/1001.4524v3. 

[47] P. Adamson, C. Andreopoulos, D. J. Auty, D. S. Ayres, C. Back- 
house, G. Barr, R. H. Bernstein, M. Betancourt, P. Bhattarai, 
M. Bishai, et al. (The MINOS Collaboration), Phys. Rev. D 82, 
051102(2010). 



